Pulse offset ultrasonic imaging

ABSTRACT

A diagnostic ultrasound image of a region of interest (ROI) of a body is formed by transmitting into the ROI at least a first and a second ultrasound pulse, in which the second pulse is phase-shifted relative to the first pulse by an amount other than 0 or a multiple of 180 degrees. Discretized receive signals from the pulses are interleaved to form a resultant operating signal that is detected and beamformed as the operating signal.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims priority of Chinese patent application no. 200910133494.X, filed 25 May 2009.

TECHNICAL FIELD

The invention relates to diagnostic medical ultrasonic imaging.

BACKGROUND ART

Ultrasonic imaging is a well-known diagnostic medical technique that can generate a visible image of a region of interest (ROI) by transmitting ultrasound into a patient's body, and then receiving and processing reflected ultrasonic echoes, whose characteristics carry information about the insonified body tissues.

In a conventional ultrasonic imaging system, a transmit circuit energizes an array of piezoelectric elements to generate ultrasonic transmit pulses along focused beams into the ROI of the body. Many body tissues of interest are typically nonlinear media such that non-linear distortion is created as the ultrasonic pulses propagate through them. The nature of this distortion indicates characteristics of the tissues that can be extracted for display.

There is a continual effort to improve the quality of these diagnostic images, with “quality” being defined in several different ways. As with many other forms of signal processing, the effort often involves trying to get as much information as possible from the received signals while rejecting as much unwanted signal energy (“noise”) as possible.

One known and common method is often referred to as “harmonic imaging.” As the name implies, instead of processing only the echo signal at the fundamental transmit frequency, this technique involves isolating and extracting information from one or more harmonics of this fundamental frequency either in addition to or instead of at the fundamental frequency alone. Such a technique is disclosed in:

-   -   U.S. Pat. No. 6,206,833;     -   U.S. Pat. No. 7,004,905;     -   U.S. Pat. No. 7,104,956; and     -   U.S. Pat. No. 7,513,870.

Echoes reflected from the ROI will typically include a linear component, which is associated with the ultrasonic pulses, and a nonlinear component, associated with the characteristic of tissues and structures of the body. The nonlinear component is extracted from the echoes, and the visible image is generated by processing and analyzing the nonlinear component.

One other known method for improving imaging quality is often referred to as “phase inversion” or “pulse inversion.” This technique exploits the linearity of much of the medium through which the ultrasound pressure waves propagate through the ROI by transmitting and receiving a first ultrasound pulse, and then soon thereafter transmitting and receiving a second ultrasound pulse with a polarity opposite that of the first. Stated differently but equivalently, the two paired pulses differ in phase by 180 degrees. The combined (typically, summed) return signals will therefore “cancel” each other out (sum to zero) with respect to all linear propagation. The remaining signal is then either processed as is, or a harmonic of it is processed to create the displayed image (thereby using both the “phase inversion” and “harmonic imaging” techniques). U.S. Pat. Nos. 5,632,277 and 5,951,478 disclose examples of systems that use the “phase inversion” technique.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram of an offset pulse ultrasonic imaging system in accordance with an embodiment;

FIG. 2 is a flowchart of an offset pulse ultrasonic imaging method in accordance with an embodiment;

FIG. 3 is a block diagram of an offset pulse ultrasonic imaging system in accordance with another embodiment;

FIG. 4 is a flowchart of an offset pulse ultrasonic imaging method in accordance with another embodiment;

FIG. 5A-5D are graphs of waveforms for ultrasonic pulses and echoes in accordance with another embodiment;

FIG. 6A-6D are graphs of echo waveforms at sample rates of 40 MHz, shown in a conventional ultrasonic imaging system and an ultrasonic imaging system in accordance with an embodiment respectively; and

FIG. 7A-7D are graphs of echo waveforms at sample rates of 30 MHz, shown in a conventional ultrasonic imaging system and an ultrasonic imaging system in accordance with an embodiment respectively.

FIG. 8 illustrates a two-pulse embodiment and how a second pulse is offset from a first.

FIGS. 9A and 9B illustrate how the pulses shown in FIG. 8 may be discretized.

FIG. 10 illustrates echo pulse interleaving.

FIG. 11 illustrates how pulse offsetting and interleaving affect the spectrum of receive signals.

DETAILED DESCRIPTION

The following detailed description refers to the accompanying drawings. The same reference numbers may be used in different drawings to identify the same or similar elements. In the following description specific details may be set forth such as particular structures, techniques, etc. in order to provide a thorough understanding of the various aspects of the claimed invention. However, such details are provided for purposes of explanation and should not be viewed as limiting with respect to the claimed invention. With benefit of the present disclosure it will be apparent to those skilled in the art that the various aspects of the invention claimed may be practiced in other examples that depart from these specific details. Moreover, in certain instances, descriptions of well known devices, circuits, and methods are omitted so as not to obscure the description of the present invention with unnecessary detail.

FIG. 1 shows a block diagram of an offset pulse ultrasonic imaging system 100 in accordance with an embodiment. As shown in FIG. 1, a transmit circuit 10 generates a plurality of ultrasonic pulse signals for each scan line in a scan plane. The plurality of ultrasonic pulse signals is converted into a plurality of ultrasonic pulses by a probe/transducer 20, which includes an array of piezoelectric elements. The structure and operation of an ultrasound probe/transducer are well understood and can be summarized by saying that it functions as a phased array of piezoelectric transmit/receive elements that generate and respond to pressure waves at ultrasonic frequencies. The ultrasonic pulses are offset relative to each other, which will be explained in detail later. The converted ultrasonic pulses are transmitted into a region of interest (ROI) 30 of a body. Since tissues of the body are nonlinear media, the echoes, reflected from the ROI 30, include a linear component associated with the ultrasonic pulses, and a nonlinear component associated with the characteristics of the various body tissues and other structures.

The ultrasound echoes reflected back from these body tissues and structures are then converted into echo signals by the probe 20. Conventional sampling, analog-to-digital and other circuitry is then used to convert the electrical signals created by the piezoelectric elements in the probe into a series of digital equivalent values, which are then supplied to an interleaver 40. As is explained in greater detail below, the interleaver 40 interleaves the values into a sequence of samples to form an operating signal. The operating signal is supplied to a detector 50. An explanation of the operating signal will be given in detail later. The nonlinear component of the is extracted from the operating signal in the detector 50 and is then beamformed in a beamformer 60 to generate a scan line, a collection of which is then converted in any know manner into a form suitable for display or recording or other processing by a user.

FIG. 2 is a flowchart of an offset pulse ultrasonic imaging method in accordance with an embodiment. In one embodiment, the method can be performed by the offset pulse ultrasonic imaging system 100 in FIG. 1.

First, at step 10, a transmit circuit 10 generates a plurality K of ultrasonic pulse signals for each scan line in a scan plane. That is, K ultrasonic pulse signals are emitted by the transmit circuit 10. For the K ultrasonic pulse signals, each ultrasonic pulse signal is time-shifted relative to an immediately preceding ultrasonic pulse signal by a predetermined amount. In one embodiment, the time-shift is τ=1/(K·f_(s)) where f_(s) is an initial sample rate for each scan line.

Each of the K ultrasonic pulse signals may optionally be weighted by a corresponding weight coefficient

${\cos \left( {\frac{{\left( {k - 1} \right) \cdot 2}\pi}{K} \cdot \alpha} \right)},$

where K=1, 2 . . . K, and α is positive integer. If the K^(th) ultrasonic pulse signal is represented by P_(K), ultrasonic pulse signals P₁, P₂, . . . , P_(K) can be represented as:

P₁ = A(t)sin (ω t + Φ_(p)) $P_{2} = {{A\left( {t + \tau} \right)}{\cos \left( {\frac{2\pi}{K} \cdot \alpha} \right)}{\sin \left( {{\omega \left( {t + \tau} \right)} + \Phi_{p}} \right)}}$ ⋯ $P_{K} = {{A\left( {t + {\left( {k - 1} \right)\tau}} \right)}{\cos \left( {\frac{{\left( {k - 1} \right) \cdot 2}\pi}{K} \cdot \alpha} \right)}{\sin \left( {{\omega \left( {t + {\left( {k - 1} \right)\tau}} \right)} + \Phi_{p}} \right)}}$

The absolute ratio of two weight coefficients represents the amplitude ratio of the two corresponding ultrasonic pulse signals. In the condition where the weight coefficient is positive, the waveform polarity of the corresponding ultrasonic pulse signal is positive; on the other hand, if the weight coefficient is negative, the polarity is negative.

At step 20, the K ultrasonic pulse signals are converted into K corresponding ultrasonic pulses by the probe 20, wherein each ultrasonic pulse is time-shifted relative to an immediately preceding ultrasonic pulse by the predetermined amount T. The K ultrasonic pulses are transmitted into the ROI 30 of the body. Because each of the K ultrasonic pulses is time-shifted relative to an immediately preceding ultrasonic pulse by T, the K ultrasonic pulses are relatively offset pulses emitted to the ROI 30.

Since many body tissues of interest comprise acoustically nonlinear media, the echoes, reflected from the ROI 30, include a linear component associated with the ultrasonic pulses and a nonlinear component associated with characteristics of tissues of the body. At step 30, the probe 20 receives K echoes and converts the K echoes into K echo signals. If the K^(th) echo signal is represented by E_(K), echo signals E₁, E₂, . . . , E_(K) (K=1, 2 . . . K) can be represented as:

     E₁ = A(t)sin (ω t + Φ_(e)) + H(t) $\mspace{79mu} {E_{2} = {{{A\left( {t + \tau} \right)}{\cos \left( {\frac{2\pi}{K} \cdot \alpha} \right)}{\sin \left( {{\omega \left( {t + \tau} \right)} + \Phi_{e}} \right)}} + {{\cos^{2}\left( {\frac{2\pi}{K} \cdot \alpha} \right)}{H\left( {t + \tau} \right)}}}}$ ⋯ $E_{K} = {{{A\left( {t + {\left( {k - 1} \right)\tau}} \right)}{\cos \left( {\frac{{\left( {k - 1} \right) \cdot 2}\pi}{K} \cdot \alpha} \right)}{\sin \left( {{\omega \left( {t + {\left( {k - 1} \right)\tau}} \right)} + \Phi_{e}} \right)}} + {{\cos^{2}\left( {\frac{{\left( {k - 1} \right) \cdot 2}\pi}{K} \cdot \alpha} \right)}{H\left( {t + {\left( {k - 1} \right)\tau}} \right)}}}$

where H(t) represents the nonlinear component which carries information on characteristic of tissues of the body.

In this embodiment, K ultrasonic pulse signals are emitted by the transmit circuit 10 for each scan line in the scan plane associated with the initial sample rate f_(s). The interleaver 40 thus may receive K echo signals converted by the probe 20 in the period of τ=1/f_(s). That is, in each sample interval of 1/f_(s), K echo signals (E₁, E₂ . . . E_(K)) can be received by the interleaver 40. These echo signals (E₁, E₂ . . . E_(K)) make up K samples in the sample interval of 1/f_(s).

Where the initial sample interval is 1/f_(s), the number of initial samples is N for each scan line in the scan plane. As mentioned above, there are K signals in the sample interval of 1/f_(s). The sampled values of the K signals in each sample interval of 1/f_(s) may then be interleaved in a sequence of the initial samples 1, 2 . . . N, to form the operating signal for generating the scan line, the number of current samples is N·K. If the operating signal is represented by W(t), W(t) may be made up of echo signals

E₁₁ E₁₂ . . . E_(1K) E₂₁ E₂₂ . . . E_(2K) . . . E_(i1) E_(i2) . . . E_(ij) . . . E_(iK) . . . E_(N1) E_(N2) . . . E_(Nj) . . . . E_(NK), i=1, 2 . . . N, j=1, 2 . . . K.

${W(t)} = {{{{A(t)}{\sin \left( {{\omega \; t} + \Phi_{e}} \right)}{\cos \left( {2{\pi\alpha}\; {f_{s} \cdot t}} \right)}} + {{H(t)}{\cos^{2}\left( {2{\pi\alpha}\; {f_{s} \cdot t}} \right)}}} = {{{A(t)}{\sin \left( {{\omega \; t} + \Phi_{e}} \right)}{\cos \left( {2{\pi\alpha}\; {f_{s} \cdot t}} \right)}} + {\frac{1}{2}{H(t)}} + {\frac{1}{2}{H(t)}{\cos \left( {4{\pi\alpha}\; {f_{s} \cdot t}} \right)}}}}$      Where  t = m τ, m = 0, 1, 2  …  (N ⋅ K − 1), τ = 1/(Kf_(s))

The K echo signals are interleaved in the sequence of the initial samples (1, 2 . . . N) to form the operating signal W(t) in the interleaver 40. The operating signal W(t) is supplied to a detector 50 at step 40. Then, at step 50, the nonlinear component is extracted from the operating signal W(t) in the detector 50, and at step 60, the nonlinear component is beamformed in a beamformer 60 to generate the scan line. Since the nonlinear component carries information about characteristics of tissues of the body, a visible image of the ROI 30 is generated by analyzing the nonlinear component using known scan conversion and display algorithms.

It can be concluded from the expression of W(t): (1) the sample interval of the operating signal is T, so the effective sample rate of the operating signal is K·f_(s) which is K times the initial sample rate f_(s); (2) the linear component of the operating signal is modulated to α times f_(s), i.e. the linear component is moved to both sides of a·f_(s) in the frequency domain; whereas, for the nonlinear component, half is not modulated and the other half is modulated to 2αf_(s). That is, in the frequency domain, half of the nonlinear component is maintained at its original frequency and the other half is moved to both sides at 2·α·f_(s).

Detection of the nonlinear component at either side of 2·α·f_(s) may avoid frequency overlap, or at least cause less overlap than in known methods. It is preferable to select α and K by considering the following equation. In accordance with Nyquist theorem:

${{KF}_{s} \geq {2 \cdot \left( \frac{4{\pi\alpha}\; f_{s}}{2\pi} \right)}},{{i.e.\mspace{11mu} \frac{K}{\alpha}} \geq 4}$

In this embodiment, each of the K ultrasonic pulse signals may be weighted by a corresponding weight coefficient, which causes the linear component to move to either side of αf_(s) in the frequency spectrum; as for the nonlinear component, half is maintained at the original fundamental frequency and the other half is moved to either side of 2αf_(s). Accordingly, when the operating signal is analyzed in the detector 50, the nonlinear component is extracted efficiently because the linear component and the nonlinear component are separated in the frequency domain.

In another embodiment, the linear component may be overlapped by part of the nonlinear component (e.g. the nonlinear component moved to 2·α·f_(s)) in the frequency domain while the nonlinear component is maintained at the original fundamental frequency (i.e. non-modulated nonlinear component) instead of the nonlinear component being moved to both sides of 2αf_(s) (i.e. modulated nonlinear component). The relevant information-bearing nonlinear signal component, can be extracted using known filtering methods provided that the non-modulated nonlinear component is sufficiently separated from the linear component and the non-modulated nonlinear component.

FIG. 3 shows a block diagram of an offset pulse ultrasonic imaging system 200 in accordance with another embodiment. As compared with the offset pulse ultrasonic imaging system 100, the system 200 further comprises a selector 5. The selector 5 may select suitable K and a for enabling the non-modulated nonlinear component to be separated from the linear component in the frequency domain and not be interfered with due to frequency overlap caused by the linear component and the modulated nonlinear component.

As shown in FIG. 3, a transmit circuit 10′, based on the selected K and a, generates a plurality of ultrasonic pulse signals for each scan line in a scan plane. The probe 20, interleaver 40, detector 50 and beamformer 60, can be implemented using the same or similar devices and programming as shown in FIG. 1.

FIG. 4 is a flowchart of an offset pulse ultrasonic imaging method in accordance with an embodiment. In this embodiment, the method can be performed by the offset pulse ultrasonic imaging system 200 shown in FIG. 3. First, at step 5, a selector 5 selects suitable K and a values. The selected K and a may enable the non-modulated nonlinear component to be separated from the linear component, thereby reducing interference due to frequency overlap of the linear and nonlinear components in frequency domain. At step 10′, a transmit circuit 10′, based on the selected K and a, generates a plurality of ultrasonic pulse signals, for example K, for each scan line in a scan plane. The following steps 20-60 may be the same or similar as those shown in FIG. 2.

The system and the method as shown in FIG. 3 and FIG. 4 will be illustrated further by providing a group of particular values of K and a. In one embodiment, K=2 and α=1. In accordance with the expression for P_(K) above, the ultrasonic pulse signals P₁ and P₂ can be represented as:

     P₁ = A(t)sin (ω t + Φ_(p)) $P_{2} = {{{A\left( {t + \tau} \right)}{\cos \left( {\frac{2\pi}{K} \cdot \alpha} \right)}{\sin \left( {{\omega \left( {t + \tau} \right)} + \Phi_{p}} \right)}} = {{- {A\left( {t + \tau} \right)}}{\sin \left( {{\omega \left( {t + \tau} \right)} + \Phi_{p}} \right)}}}$

In this example, the polarities of waveforms of ultrasonic pulse signals P₁ and P₂ are positive and negative respectively, and the amplitude ratio of P₁ and P₂ is 1:1. The ultrasonic pulse signal P₂ is time-shifted relative to the ultrasonic pulse signal P₁ by a predetermined amount σ=1/(2f_(s)). The waveforms for the ultrasonic pulses P₁ and P₂ are illustrated in FIG. 5A.

The ultrasonic pulses P₁ and P₂ are transmitted into the ROI 30 and echo signals are reflected back to the transducer as usual. In accordance with the expression of echo signal E_(K) above, the echo signals E₁ and E₂ can be represented as:

     E₁ = A(t)sin (ω t + Φ_(e)) + H(t) $E_{2} = {{{{A\left( {t + \tau} \right)}{\cos \left( {\frac{2\pi}{K} \cdot \alpha} \right)}{\sin \left( {{\omega \left( {t + \tau} \right)} + \Phi_{e}} \right)}} + {{\cos^{2}\left( {\frac{2\pi}{K} \cdot \alpha} \right)}{H\left( {t + \tau} \right)}}} = {{{- {A\left( {t + \tau} \right)}}{\sin \left( {{\omega \left( {t + \tau} \right)} + \Phi_{e}} \right)}} + {H\left( {t + \tau} \right)}}}$

The waveforms for the echo signals E₁ and E₂ are illustrated in FIG. 5B, where the echo signal E₁ is illustrated with small circles and the echo signal E₂ is illustrated as small squares. The two echo signals E₁ and E₂ in each sample interval of 1/f_(s) are interleaved to form an operating signal W(t). In accordance with the expression of the operating signal W(t) as above mentioned, W(t) can be represented as:

W(t)=A(t)sin(ωt+Φ _(e))cos(2πf _(s) ·t)+½H(t)+½H(t)cos(4πf _(s) ·t)

In this embodiment, the sample rate of the nonlinear component is 2f_(s) (K=2), so the operating signal W(t) can be further expressed as:

W(t)=A(t)sin(ωt+Φ _(e))cos(2πf _(s) ·t)+H(t)

It can be concluded from the expression of W(t): (1) the sample rate of the operating signal is two times of the initial sample rate f_(s), i.e. 2·f_(s), as shown in FIG. 5C; (2) the linear component is moved to either side of f_(s) in frequency domain, the two central frequencies are located at (f_(s)−ω/2π) and at (fs+ω/2π) respectively; and the nonlinear component is maintained at the fundamental frequency. Since the linear component and the nonlinear component are separated in the frequency domain, and the sample rate of the nonlinear component is two times the initial sample rate f_(s), the linear component can be efficiently filtered using known filtering methods, for example, and the nonlinear component H(t) can be extracted efficiently as shown in FIG. 5D.

Further, in one illustrative embodiment, the initial sample rate f_(s) is 40 MHz, the central frequency of the linear component is 6 MHz and the central frequency of the nonlinear component is 12 MHz. The waveforms of the operating signal W(t) are illustrated in FIG. 6C and FIG. 6D; whereas the waveforms of echoes in a conventional ultrasonic imaging system are shown in FIG. 6A and FIG. 6B. There will be significant frequency overlap in conventional ultrasonic imaging systems; however, according to the embodiment of the invention (K=2, α=1), not only the effective sample rate is increased to two times the initial sample rate f_(s) (i.e. 80 MHz), but the linear component is much more completely separated from the nonlinear component in the frequency domain.

According to embodiments of the present invention, the sample rate of the nonlinear component is increased by virtue of the parameter K. Thus, almost all frequencies in the frequency band of the nonlinear component can be more efficiently detected even if the initial sample rate f_(s) is not high.

In another embodiment, the initial sample rate f_(s) is 30 MHz. The waveforms of the operating signal W(t) are shown in FIG. 7C and FIG. 7D; the waveforms of echoes in a conventional ultrasonic imaging system are shown in FIG. 7A and FIG. 7B. As can be seen, because of significant frequency overlap in the conventional ultrasonic imaging system, the information at the frequencies in the frequency band of the nonlinear component cannot be detected completely; however, according to the embodiment of the present invention, not only the linear component can be separated from the nonlinear component in the frequency domain, but the sample rate can be increased to two times of the initial sample rate f_(s) (in this example, 60 MHz), which enables almost all signal at frequencies in the frequency band of the nonlinear component to be much easier to isolate and detect.

Other groups of particular values of K and a can be selected. With the selected K and a, the effective sample rate of the operating signal is increased to K times the initial sample rate f_(s), and the linear component is modulated to a location as far as possible from the nonlinear component in the frequency domain.

According to another illustrative embodiment, K=6, α=2, such that six ultrasonic pulse signals P₁, P₂ . . . P₆ are generated for each scan line in the scan plane. The corresponding weight coefficients of the six ultrasonic pulse signals may then be 1, −½, —½, 1, −½, —½. The polarities of the waveforms for the six ultrasonic pulse signals are thus positive, negative, negative, positive, negative, negative, and the amplitude ratios of the six ultrasonic pulse signals are 1:0.5:0.5:1:0.5:0.5. For the six ultrasonic pulse signals, an ultrasonic pulse signal is time-shifted relative to an immediately preceding ultrasonic pulse signal by a predetermined amount τ=1/(2f_(s)).

The ultrasonic pulses P₁, P₂ . . . P₆ are transmitted into the ROI 30, and the echoes reflected from the region 30 and sensed by the transducer as usual. The six echo signals E₁, E₂ . . . E₆ are then interleaved in a sequence of initial samples to form an operating signal W(t). In accordance with the expression of the operating signal W(t) as above, W(t) can be represented:

W(t)=A(t)sin(ωt+Φ _(e))cos(4πf _(s) ·t)+½H(t)+½H(t)cos(8πf _(s) ·t)

It can be concluded from the expression of W(t): the linear component is moved to either side of 2·f_(s) in frequency domain; the two central frequencies are located at (2f_(s)−ω/2π) and (2f_(s)+ω/2π) respectively; for the nonlinear component, half of which is maintained at the original fundamental frequency and the other half of which is moved to either side of 4·f_(s); and the two central frequencies are located at (4f_(s)−ω/π) and (4f_(s)+ω/π). In this embodiment, the effective sampling rate of the operating signal is increased to six times the initial sampling rate, that is, to 6·f_(s) (K=6). Due to the interference of frequency overlap, within an active range of 0˜3 f_(s) for frequency detection, the linear component, which is moved to both sides of 2·f_(s), is overlapped or partly overlapped with the nonlinear component, which is moved to both sides of 4f_(s) (−4f_(s)+6f_(s)=2f_(s)) at 2·f_(s) in the frequency domain. However, the other half of the nonlinear component, which is maintained at original fundamental frequency, is not significantly affected by the frequency overlap. Since half of the nonlinear component, which is maintained at original frequency, is separated from the linear component, which is moved to either side of 2·f_(s) in the frequency domain, and the half of the nonlinear component, which is maintained at original frequency, is not interfered with significantly by the frequency overlap, the linear component can be filtered out using conventional filtering methods and the nonlinear component H(t) can be extracted from the operating signal W(t).

In an ultrasonic imaging system according to embodiments of the invention, K ultrasonic pulse signals are emitted into a region of interest of a body for each scan line in a scan plane. Accordingly, K echoes are reflected by the region of interest. An operating signal is formed by interleaving K echoes signals. When the operating signal is detected to extract the information from the nonlinear component of the echoes, the frequency of the linear component of the echoes is modulated away so as to enable more efficient separation from the nonlinear component in frequency domain.

The offset pulse ultrasonic imaging system and method thereof should not be limited to embodiments mentioned above. It will be apparent to those skilled in the art that the various aspects of the invention claimed may be practiced in other examples that depart from these specific details.

In an alternative embodiment of the ultrasonic imaging method, a plurality of groups of ultrasonic pulses are transmitted into a region of interest of a body for each scan line in a scan plane, wherein each group of ultrasonic pulses comprises K ultrasonic pulses. For example, a transmit circuit may generate three groups of ultrasonic pulses for each scan line in a scan plane, each group of ultrasonic pulses comprising two ultrasonic pulses P₁ and P₂, wherein P₂ is time-shifted relative to P₁ by a predetermined amount, and the two ultrasonic pulses P₁ and P₂ are optionally weighted by corresponding weight coefficients. The predetermined amount and the corresponding weight coefficients may adopt parameters disclosed in above embodiments. The three groups of ultrasonic pulses (P₁, P₂, P₁, P₂, P₁, P₂) are emitted into the region of interest. The echoes reflected from the region of interest are received by the probe/transducer as usual. The received echoes accordingly include six echo signals (E₁, E₂, E₁, E₂, E₁, E₂). Each group of echo signals (E₁, E₂ pair) is interleaved in a sequence of initial samples to form an operating signal associated with the group of echo signals. Therefore, three groups of echo signals (E₁, E₂ pairs) are interleaved to form three operating signals W₁, W₂, W₃. The nonlinear component associated with each of the groups of echo signals is then extracted from the corresponding operating signal; that is, nonlinear components h₁, h₂ and h₃ are extracted from W₁, W₂, W₃ respectively. In a module such as the detector in above embodiments, the extracted nonlinear components h₁, h₂ and h₃ can be combined, for example, as a resultant of the extracted nonlinear components, by summing the nonlinear components, i.e. H=h₁+h₂+h₃. The SNR (signal-to-noise ratio) of H may thereby be improved in comparison with the SNR of each of linear components h₁, h₂ and h₃.

In addition, the selector described above may be configured to allow for user adjustment. The selector may select suitable K and a values based on the active range of frequency detection. With the selected K and a, a transmit circuit generates K ultrasonic pulse signals for each scan line in a scan plane. In an alternative embodiment, a particular transmit circuit may generate a predetermined plurality of ultrasonic pulse signals for each scan line in a scan plane.

Further, as in other embodiments, the amplitude of each of the K ultrasonic pulse signals may be weighted by a corresponding weight coefficient. In an alternative embodiment, K ultrasonic pulse signals are emitted to a region of interest of a body for each scan line in a scan plane, wherein the amplitude of each of the K ultrasonic pulse signals is not weighted. K echoes are reflected by the region of interest accordingly. An operating signal is formed by interleaving K echoes signals. When the operating signal is detected to extract nonlinear component of the echoes, the effective sampling rate of the operating signal is increased to K times the initial sample rate f_(s). Since the linear component is not modulated, the linear component may be overlapped with the nonlinear component in the frequency domain. A narrow-band filter that passes narrowband ultrasonic pulse signals may then be used to extract the nonlinear component in this embodiment.

Also, an ultrasonic pulse signal is time-shifted relative to an immediately preceding ultrasonic pulse signal by a predetermined amount. In some of the above-mentioned embodiments, the time-shift τ=1/(K·f_(s)) wherein f_(s) is the initial sample rate for each scan line. In an alternative embodiment, the time-shift τ may be set another predetermined value, for example 1/(2·K·f_(s)). Correspondingly, some modifications might be made to the ultrasonic imaging system. For example, a synchronizer may be added into the system as needed to keep the ultrasonic pulse signals emitted by the transmit circuit and the received echoes synchronized so as to form an operating signal.

FIGS. 8-11 illustrate in a different and perhaps more visible way the two-pulse embodiment of the invention described above. The illustrated example also shows the principle of the invention in terms of phase shift instead of time shift, although these two representations are of course equivalent mathematically. In the illustrated example, two transmit pulses P1(t) and P2(t) are shown in greatly simplified form for the sake of clarity; all skilled in medical ultrasound engineering realize that actual transmit pulses from a modern ultrasound imaging transducer—especially with the usual multi-element phased piezoelectric array—may have a much more complicated waveform.

The main point to be illustrated in FIG. 8 is that the second pulse in the series is generated so as to be offset in phase from P1 by Δθ, which is neither 0 degrees nor 180 degrees. In other words, in this invention, the second pulse P2 is not just a “repeat” (0° phase shift or a multiple of 360°) of the first pulse P1, nor is it the “inverse” (180° phase shift or a multiple thereof) of P1, but is rather offset from P1 by an amount that is other than 0° or 180°.

Sampling of analog signals to generate corresponding numerical representations is of course a well-known feature of almost all digital signal processing. In FIGS. 9A and 9B, receive/echo signals E1, E2 corresponding to transmit pulses P1, P2, respectively, are illustrated in discretized form, such that E1=E1(i), i=0, 1, . . . and E2=E2(j), j=0, 1, . . . .

FIG. 10 illustrates graphically how E1 and E2 may be interleaved such that the values E2(j) are aligned “in between” temporally adjacent values of E1(i). If P2 is sampled at the same frequency as P1, then each E2(j) value can be interleaved temporally in the middle of the adjacent E1(i) values. The combined, that is, interleaved resultant operating signal W will thereby have twice as many values per cycle as either E1 or E2 alone; in other words, the effective sampling rate of the interleaved resultant signal is twice the actual sampling rate used to discretize each of the return signals from P1 and P2. The combined, interleaved operating signal W can then be processed as if it were a single discretized receive signal

Those skilled in the art will readily see how the two-pulse embodiment illustrated in FIGS. 8-10 can be extended to even more pulses with suitable relative offsets. In other words, although FIGS. 8-10 illustrate an example with K=2, the principle can be extended to embodiments in which K>2.

FIG. 11 illustrates the effect of transmit pulse offsetting and receive pulse interleaving. Assume first that there were no such offsetting and interleaving. The spectral component at the fundamental frequency f₀ (shown with the dashed line) and the spectral component at the first harmonic (centered on 2f₀) will overlap (shown with hatched shading); moreover, even a narrow-band filter will typically also pass energy from the fundamental component in addition to the harmonic information, noise, and non-linear information N. One prior art way to try to overcome this inherent problem is to then try to extract even higher-order harmonic information, but as is well known, the amplitude of such higher-order harmonics becomes lower and lower, thereby making it more and more difficult to maintain a good signal-to-noise ratio.

The effect of transmit pulse offsetting and return signal interleaving and compounding is that the spectral component at the fundamental frequency f₀ is shifted both upward and downward relative to the origin, with the positive shift (illustrated by the arrow arc from the dashed f₀ “curve” to the solid f_(s)+f₀ curve. One advantage of this is that a band-pass filter narrowly tuned to the harmonic 2f₀ is not necessary; rather, even a broader-spectrum filter F can suffice; alternatively, the signal from a band-pass filter will not be as degraded by component overlap.

The amount of offset, that is, Δθ may be adjusted according to the needs of any given implementation of the invention and could even be made user-adjustable during scans to find a “best” image result. In general, as AO increases, the degree of separation, that is, fundamental component up-shift, will increase, but accuracy will decrease. For K=2, a suitable value for Δθ will typically be Δθ=(f_(s)/2f₀)( 1/360). 

What is claimed is:
 1. A method for creating an image of a region of interest of a body comprising: transmitting a first ultrasonic pulse into the region of interest; sensing first receive signals corresponding to ultrasonic energy from the first ultrasonic pulse reflected from within the region of interest; transmitting a second ultrasonic pulse into the region of interest; sensing second receive signals corresponding to ultrasonic energy from the second ultrasonic pulse reflected from within the region of interest; interleaving discretized representations of the first and second receive signals to form a resultant operating signal; and creating the image of the region of interest as a function of the resultant operating signal; wherein the second ultrasonic pulse is phase-shifted relative to the first ultrasonic pulse by an amount other than zero degrees or a multiple of 180 degrees.
 2. The method of claim 1, further comprising: extracting a nonlinear component from the operating signal; and beamforming the nonlinear component to generate a scan line, whereby a series of scan lines form the image.
 3. The method of claim 1, in which K receive signals, corresponding to echo signals from respective, relatively phase-offset transmitted ultrasonic pulses, are discretized and interleaved to form the resultant operating signal.
 4. The method of claim 3, in which K is greater than two.
 5. The method of claim 3, in which: K is two or greater; the phase shift is chosen such that the second ultrasonic pulse is transmitted with a time delay of 1/K·f_(s); and the receive signals are discretized with an initial sampling rate of f_(s).
 6. The method of claim 3, wherein amplitudes of the K ultrasonic pulses are weighted by corresponding weight coefficients.
 7. The method of claim 6, wherein the weight coefficient for the amplitude of a k^(th) ultrasonic pulse is ${\cos \left( {\frac{2{\pi \left( {k - 1} \right)}}{K} \cdot \alpha} \right)},{k = 1},{2\mspace{14mu} \ldots \mspace{14mu} K},$ α is positive integer and the center frequency of the linear component is moved to αf_(s).
 8. The method of claim 7, wherein the (K, α) is selected from any one of groups of (2, 1) and (6, 2).
 9. The method of claim 2, further comprising frequency filtering the operating signal using a filter than passes not only a frequency that is a harmonic of a fundamental signal frequency, but also energy at the fundamental signal frequency itself.
 10. A system for creating an image of a region of interest of a body comprising: an ultrasound transducer transmitting a first ultrasonic pulse and at least a second ultrasonic pulse into the region of interest, the ultrasound transducer sensing first and second receive signals corresponding to ultrasonic energy from the first and ultrasonic pulses, respectively, from within the region of interest; transmit control circuitry for causing the ultrasound transducer to generate the second ultrasonic pulse phase-shifted relative to the first ultrasonic pulse by an amount other than zero degrees or a multiple of 180 degrees; an interleaver for interleaving discretized representations of the first and second receive signals to form a resultant operating signal; and detection and beamforming circuitry for creating the image as a function of the resultant operating signal.
 11. The system of claim 10, in which the detection and beamforming circuitry comprises: a detector to extract a nonlinear component from the operating signal for each group of ultrasonic pulses and to generate a resultant of the nonlinear component for each group of transmitted ultrasonic pulses; and a beamformer to beamform the resultant of the nonlinear component for each group of ultrasonic pulses to generate a scan line, whereby a series of scan lines form the image.
 12. The system of claim 10, in which K receive signals, corresponding to echo signals from respective, relatively phase-offset transmitted ultrasonic pulses, are discretized and interleaved to form the resultant operating signal.
 13. The system of claim 12, in which K is greater than two.
 14. The system of claim 12, in which: K is two or greater; the phase shift is chosen such that the second ultrasonic pulse is transmitted with a time delay of 1/K·f_(s); and the receive signals are discretized with an initial sampling rate of f_(s).
 15. The system of claim 14, wherein amplitudes of the K ultrasonic pulses are weighted by corresponding weight coefficients.
 16. The system of claim 15, wherein the weight coefficient for the amplitude of a k^(th) ultrasonic pulse is ${\cos \left( {\frac{2{\pi \left( {k - 1} \right)}}{K} \cdot \alpha} \right)},{k = 1},{2\mspace{14mu} \ldots \mspace{14mu} K},$ α is positive integer and the center frequency of the linear component is moved to αf_(s). 